ERROR ESTIMATES FOR A NEUMANN PROBLEM IN HIGHLY 
OSCILLATING THIN DOMAINS 
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Abstract. In this work we analyze the convergence of solutions of the Poisson equation 
with Neumann boundary conditions in a two-dimensional thin domain with highly oscillatory 
behavior. We consider the case where the height of the domain, amplitude and period of 
the oscillations are all of the same order, and given by a small parameter e > 0. Using an 
appropriate corrector approach, we show strong convergence and give error estimates when 
we replace the original solutions by the first-order expansion through the Multiple- Scale 
Method. 



Several important problems arising in physics and engineering lead to consider boundary- 
value problems in domains with oscillating boundaries, such as flows over rough walls, elec- 
tromagnetic waves in a region containing a rough interface [2D], or heat transmission in 
winglets [2]. 

In the present paper we consider the Poisson equation with Neumann boundary conditions, 
arising from the study of stead state of reaction-diffusion flow over a thin bar with a rough 
boundary [1]. The aim of our study is to derive precise estimates for the effective behavior 
of the solutions of such problem. 

In order to setup the problem, let g G C^(]R, M) be a L-periodic positive function. Given 
a small parameter e > 0, let us consider the family of two-dimensional domains 



where N"^ denotes the unit outward normal vector field to dR' and is a non-homogeneous 
term in L'^{R'^) uniformly bounded. 
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1. Introduction 



i?' = {{xi,X2) G I < xi < 1, < a;2 < eg{xi/e)}. 
We restrict our attention to the solutions w' of the family of elliptic equations 

f -Aw' + w' = r in R' 





(1.2) 
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Figure 1. The thin domain i?^ 



Many authors have devoted efforts in the investigation of the asymptotic behavior of a 
similar equation in related domains, e.g. [21 El [131 IIHI ED]- The basic assumption is that 
the scale of wavelength of the roughness is small compared with the roughness amplitude. 
Therefore it is important to emphasize that the amplitude and period of the oscillations of 
are of the same order e, which also coincides with the order of thickness of the thin domain. 
This scaling makes the problem very resonant and the determination of the limiting problem 
(e = 0) is not straightforward. It makes our analysis and results different from those papers. 



Going back to problem (1.2), in [4J the authors combine methods from homogenization 
theory, specially those related to reticulated structures, to obtain the limiting behavior of 
the family of solutions w^. We recall that the convergence obtained in ^ is with respect to 
the weak topology of the Sobolev space H^. As pointed in |3j, this convergence is the best 
one in these spaces. In general, one can not expect strong convergence in homogenization 
related problems. We refer the reader to [H HU EZj for a general introduction to the 
theory of homogenization and to [T^ for a general treatise on reticulated structures. 

In such problems, where the analytical solution demands a substantial effort or is unknown, 
it is of interest in discussing error estimates (in norm) when we replace it with numerical 
approximat ions . 

In this work we made use of an appropriate corrector approach developed by Bensoussan, 
Lions and Papanicolaou in [7] to derive a kind of strong convergence and error estimates in 
if ^-norm. This can be summarized as: 

Suppose w'' Wq weakly in H^i^R") in a sense to he defined. It is possible to find an 
explicit corrector term G H^{R'^), = o(e) in L^{R^), in order to 



for e > small enough. The reader is invited to [TO] for a classical introduction to correctors 
approach, [T71 [181 121] for further discussion on approximations of arbitrary order for related 
problems, and [HI [23 for recent works in this subject as well. 

Still related with the limiting behavior of problems in thin domains but in connection 
with continuity properties of global attractors associated to parabolic problems, we can cite 
[H [151 [2211231 [21 [26] and their references. 

Finally, we would like to observe that although we just treat the Neumann problem for 
the Poisson equation, we may also consider different conditions in the lateral boundaries of 



ERROR ESTIMATE FOR A NEUMANN PROBLEM 



3 



the thin domain R^, while preserving the Neumann type boundary condition in the upper 
and lower boundary. The limiting problem preserves the boundary condition. Note that 
Dirichlet boundary condition in the upper and lower boundaries is a trivial case since the 
trace operator is continuous on the set (0, 1) C 

This paper is organized as follows. In Section |2j we introduce the functional setting 
for the perturbed problem (1.2) and the limiting one. In Section Isl we provide formal 



computations to derive the homogenized problem by the Multiple- Scale Method. In Section 
|4| we obtain convergence in if^-norms for the first-order corrector, and in Section [sj we use 
the second-order corrector to obtain rates of convergence for the first order approximation. 
We emphasize the relevance of these rates for numerical approximations of solutions of partial 
differential equations in highly heterogeneous and porous media. 



2. Preliminaries 

We stress for the fact that R"^ varies in accordance with a positive parameter e and, when 
e goes to 0, the domains i?^ collapse themselves to the interval (0,1). Therefore, in order 
to preserve the "relative capacity" of a mensurable subset O G R'^, we rescale the Lebesgue 
measure by 1/e, dealing with the singular measure 

p,(0)=e"i|0|. 

This measure has been considered in studies involving thin domains, e.g. [T5| [22| [21], and 
allows us introduce the Lebesgue L'^{R^; Pe) and the Sobolev H^{R''; spaces. The norms 
in these spaces will be denoted by ||| • |||L2(^e) and ||| ■ respectively, being induced by 

the inner products 

{u,v)e = e^^ / uvdx, y u,v E L'^{R'') 

and 



ae{u,v) 



,-1 



e 



/ {Vu-Vv + uv}dx, Wu,v e H^{R') 
Jr" 



respectively. 



Remark 2.1. Notice that the \ \\ ■ \ \\ - norms and the usual ones in L'^{R'^) and H^{R'^) are 
equivalents and easily related by 

\\H\\l2{r^) = e"^''lM||L2(i?^), Vm e L'^{R^) 
\\\u\\\m{Rn = e'^^^lHmiR^, Vm g H\R'). 



The variational formulation of (1.2) is: find G H^{R'^) such that 

/ \Vw' ■V<^ + w'Adxidx2= I f<^dxidx2, V<^ G if^(i?'), (2.1) 

Jr^ ^ ^ Jr'^ 

which is equivalent to find w'^ G H^{R'^] p^) such that 

a,{^M) = {^J%. Vy.G//i(i?^p,). (2.2) 
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Observe that the solutions w'^ satisfy a priori estimates uniformly with respect to e. Thus, 
we can take ip = in (2.1) and (2.2) to obtain 



|vw'iiii2(^.) + iik'iiii2(H.) < \\\r\\\L^R^)\\H\\y(^R^) 



In order to capture the limiting behavior of a^{w'^, w'^) as e — )■ 0, we consider the sesquilinear 
form ao in i?^(0, 1) given by 



^u,v &H\<d,l), (2.4) 

where g = g{s) ds is the average of function g on the interval (0,L). We also will 
consider L^(0, 1) endowed with the norm induced by the inner product (■, ■)o, given by 

{u,v)(j = g / uvdx, Vm, f G L^(0, 1). (2.5) 



3. MuLTiscALE Expansion of Solutions 

In this section we proceed as in |3l El [TUl [121 We suppose for a moment /""(xi, X2) = 
/(xi) for all e > and (xi,X2) G -R^. Then we use the Multiple-Scale Method to introduce 
formally the homogenized equation and the second-order expansion to the solutions G 
if^(i?^) of problem (1.2). We seek for a formal asymptotic expansion of the form 

w {xi,X2) = Wo {^xi, —, —j + ewi {^xi, — , + e W2 [xi, — , — j + . . . (3.1) 

The variables Xi and X2 represent the "macroscopic" scale on the model, while Xi/e and X2/e 
represent the "microscopic" effect of the oscillating phenomena on the thin domain. 

The fact that degenerates to a line segment when e goes to suggests that the solutions 
w'^ will not depend on the "macroscopic" variable X2- This is taken into account in the choice 
of Wi assuming that Wi does not depend on the macroscopic variable X2. 

In order to construct the functions Wi of (3.1), motivated by the periodic nature of R", we 
consider the basic cell 

Y* = {{y, z) G I < y < L, < 2 < g{y)}. (3.2) 

We decompose dY* in the lateral part of the boundary Bq = {{0, z) & \ < z < g{0)} U 
{(L, z)) G I < 2 < g{L)}, the upper part of the boundary Bi = {{y, g{y)) G | < y < 
L} and the lower part of the boundary i?2 = {(y, 0)) G -R^ | < ?/ < L}. We also assume that 
Wi = Wi{x,y,z) is a L-periodic function in the variable y (ie. Wi{x,y + L,z) = Wi{x,y,z)) 
which is defined for all x G (0, 1) and {y, z) G F*. 

Performing the change of variables x = xi, y = — , z = — , we obtain that 

(9 = (9 + -(9 

2 ' I '1 (3-3) 

9x1X1 9xx ~l~ ~dxy ~^ ~^9yy^ 9x2X2 ~^9zz- 



Xi 


_ X2 


— , z = 




e 




9x2 = 
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Observing from the geometry of that the unit outward normal vectors to dW and to dY*, 
= [Nl, N2) and = (A''i, A^2) respectively are related by 

'Xi 

e e 



N\xi,X2)=N 



a.e. in dR^ 



(3.4) 



one can see by replacing the expansion (3.1) in the problem (1.2), after equate powers of e, 
that the differential equations for the functions wq, wi and W2 are as follows: For wq 

-Ay,2Wo(x, = in Y* 



dN 



[x,y,z) 



on 5i U B2 



Wo{x, - jz) L - periodic 
which implies that Wq has to be constant with respect to variables z, ie, 

Wo{x, y, z) = Wo{x), V {x, y, z) G (0, 1) x y*. 
Concerning Wi, we have 

-/^y^;,wi{x, y, z) = 2d^dyWo in Y* 



dwi 
'dN 



.x,y,giy)) = 
dwi 



g'iy) 



dwn 



X, y,z) = on B2 



on Bi 



dN 

wi{x, ■,z) L - periodic. 



Since dxdyWo = 0, denoting by X{y, z) the unique solution (up to a constant) of 

-A,,,X(y,^) = inF* 

on Bi 

^l + {g'{y))- 

dX 



dX g\y) 



we get 



^^-(l/,^) = onB2 
X{-, z) L - periodic, 



wi{x,y,z) = -X{y,z)^{x). 



For W2, the boundary value problem is 

' -Ay^^W2ix,y,z) = f{x) -woix) + 2dxdyWi{x,y, z) + dlwo{x) in Y* 



dw2 



[x,y,g{y)) 



g'{y) dwi 



x,y,g{y)) on Bi 



dwo 



dN 

W2{x, -, z 



X, y,z) = on B2 
L - periodic. 



(3.5) 



(3.6) 



(3.7) 



(3.8) 
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or, equivalently, upon using the expression for wi in (3.8), 



-Ay^^W2{x,y,z) = f{x) - Wo{x) + (1 - 2dyX{y, z))——{x) in Y 



dw2 
'dN 



ix,y,giy)) 



dx"^ 

X[y,g[y))^^(x) 



dx"^ 



on Bi 



(3.9) 



= on ^2 

^2(0;, - jz) L - periodic. 

Remark 3.1. These problems are second order partial differential equations in the variables 
{y, z) G y*, with x G (0, 1) playing the role of a parameter. 

Next, by Fredholm alternative we can derive the hmiting problem for wq- Indeed, if we 
take test functions (j){x,y,z) = (j){x) in (3.9), we get 







Y* 



+ 



(x) -wo(x) + {l-2dyX{y,z))-^{x) \dydz 



(3.10) 



On the other hand, by the Divergence's theorem and (3.9) 

72 



0(x; 



d Wq 
dx"^ 



x)dyX{y, z)dydz 











Jy* 






/ 

JdY* 


'p{x) 




/ 

JdY* 


'p{x) 



dx"^ 

dx'^ 
dw2 



x)Vy^-,X{y, z) • 1^ J ^ dydz 
'x)dyX{y,z)NidS 



ON 



{y, z) dS. 



Therefore replacing this expression in (3.10), we reach 



x) \ fix) - Wo{x) + (1 - dyX{y, z)) -^{x) \dydz, V0 G C^{0, 1). (3.11) 



Since X depends only on y and z, we can conclude that Wq must satisfy 

d'^w 
dx"^ 



-r-^{x) +Wq{x) = f{x), xG(0, 1) 







where 



\Y 



<(0) = w',{l] 
^\ I ~ ^^'^*^^'^)}^^^^• 



(3.12) 



(3.13) 



Remark 3.2. The second order differential equation (3.12) with the constant coefficient r is 
called the homogenized equation of problem (1.2) with the homogenized coefficient r . 
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Now we use the homogenized equation (3.12) to describe W2- We rewrite problem (3.9) as 



-divj^,^ Vy,zW2ix,y,z) 







{l-r-dyX{y,z))-^{x) 



dw2 , d'^Wo 

{x,y,9{y)) = X{y,g{y))-^{x)Ni on Bi 



dN 



dx"^ 

{x, y,z) = on B2 
W2{x, z) L - periodic. 



(3.14) 

The hnearity of (3.14) together with the fact that does not depend on the variables y 
and z suggest that we look for W2{x, y, z) of the form 



W2{x, y, z) = 0{y, z) -^i.^) for ^ ^ (0, 1) and {y, z) e Y* 



(3.15) 



where 6 is the solution of the auxiliary problem 

X{y,z) 




-divj^,^ \^y^z9{y, z) 
Vy^z9{y,z) 



1 - r - dyX{y, z) in Y* 



X{y,z) 




■ = on 5i U ^2 
9{-, z) L - periodic. 



(3.16) 



Then, we can use (3.8) and (3.15) to introduce the following asymptotic expansion for (1.2): 

'Xi X2\ d'^Wo 



w^ixi, X2) = Wo{xi) - eX 



Xi X2\ dw, 



e e 



{xi) + eH(^- 



dx"^ 



■(xi) + 



(3.17) 



dx \ e e 

This expansion will play an essential role in Section |5] below. 

Remark 3.3. According to Bensoussan, Lions and Papanicolaou in [7J, the functions X and 
9 define the first-order correctors 

'Xi X2\ dwo 



K^{xi,X2) = -eX ( — , —] -f^ixi), (xi,X2) e -R' 



and the second-order correctors 

'Xi X2\ dwo 



/i'(Xi,X2) 



-eX 



e e 



e e / dxi 

'Xi X2\ d'^Wo 



(3.18) 



^^(xO + e^^(i^,^) ^(xO, (x„X2)Gi^^ (3.19) 

Remark 3.4. The functions X and 9 are originally defined in the representative cell Y* , but 
to consider these functions in the thin domain R^, we use their periodicities at y to extend 
them to the hand 

r = {(y,z)GM'|yGM,0<z<(7(y)}, 
and we compose them with the diffeomorphisms 

T' -.R'^Y : (xi,X2) ^ (xi/e,X2/e). 
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In the analysis below, with some abuse of notation we will denote these compositions by 
X {xi/e,X2/€) and 9 {xi/e,X2/^) everywhere for {xi,X2) G -R^ With these considerations we 
can obtain some estimates on for X and 9 as well. It is easy to see that 



\^ \\l'^{r^) 



|X(xi/e,a;2/e)|^ dxidx2 



l/eL 
fc=l 

<e/L||X"^ 



\X{y,z)fdydz 



(3.20) 



2 

L2(y.)- 



Similarly we can get 







< 


e/L\\9\\l,^y.), 


dyX 


Il2(R.) 


< 


e/L\\dyX\\l,^y,^ 


d.X 


|2 

\l-^{R') 


< 


e/L\\d.X\\l,^y.) 


\\dy9 


|2 

\l'^{r^) 


< 


e/L\\dy9\\l2^Y'-)^ 


\\d,9 


|2 

\L^{R^) 


< 


e/L\\dMlHY'y 



(3.21) 



Remark 3.5. We can solve the problems (3.5), (3.6), (3.7), (3.14) and (3.16) applying the 
Lax-Milgram Theorem to the elliptic form 



ay* if A) = '^y,z<f ■'^y,z(pdydz, \f!f,(j)eH^{Y*) 
Jy* 

on the set V = Vy* /M where 

Vy, = {^e H\Y*) I (f is L- periodic with respect to y-variable}. 



Indeed, the following quantity 



defines a norm on V. 



M\v 



dydz 



1/2 



Remark 3.6. Also, we can use the elliptic form ay* to show that the homogenized coefficient 
r is positive. We will perform this here for reader's convenience. For all (p & V , we have 
that the solution X satisfies 



ay^X, 



Ni(PdS. 



Bi 



Recall that Bi is the upper part of the boundary of the basic cell. Consequently, yi — X 
satisfies 



ay4y-X,(f))= Ni^dS- Ni(PdS = 



Bi 



(3.22) 



Bi 
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for all (J) E V . Also, we have by (3.13) 

r\Y*\ = j ^^{y-X{y,z))^dydz = j V{y - X{y, z)) ■ Vy dydz 

= aY*{y-X,y). (3.23) 
Hence, due to relation (3.22) with cf) = —X and identity (3.23), we get 

r\Y*\ = ay^y - X,y) + ay^iy - X,-X) = ay^iy - X,y - X) > 0. (3.24) 

4. First-order Corrector 



As already noted, the solutions w"^ of ( 1.2 ) actually do not converge in if^-norms. However, 
if we "improve" w"^ by its first-order corrector, we are able to show the following result 

Theorem 4.1. Let w"^ be the solution of problem (1.2) with G L'^{R'^) satisfying 

|||/1l|L2(i?^) < C 

for some C > independent of e. Consider the family of functions G L'^{0, 1) defined by 

"egixi/t) 



r(xi) 



f^{Xi,X2)dX2. 



\im\\\w' - Wo - k'\\\hi(r^) = 0, 



is the unique solution of the homogenized equation (3.12) with 



(4.1) 
(4.2) 



//f w-L\0,l), then 

where k'^ is the first-order corrector of defined in Remark 3.3 and wq G H'^{0, 1) nC"'^(0, 1) 

(4.3) 



/o = T /. 
9 



Proof. By variational formulation of (1.2), we have 

a,{v,wn = {^,f%, \/veH\W). 
Thus, observing that Wq + E H^{R'^^ we obtain by symmetry of 

\\\w'' — Wq — = a^{w'' — Wq — , W^ — Wq — k") 

= a^iw" .w" — Wq — K^) — a^iwQ + k'^jW^) + a^{wQ + k^, Wq + k") 

= {w' - 2{wq + k'), /'), + a^iwQ + k\ Wq + k"). 

(4.4) 

Using the change of variables {x,y) — )■ {x,y/e) on [4, Theorem 4.3], it is easy to see that 

e~^||w'' — Wq\\l2(^ji,) — )• as e — )■ 0. 

e-S>0 



Consequently, 1 1 — Wo| | |l2 



0. Therefore 



{w^-WqJ^), ^ \\\w' - Wq\\\l2 



-)■ 0, as e 0. 



(4.5) 



* Here wq is considered as a function of xi and X2, simply with some abuse writing wo{xi,X2) 
wo{xi). 
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By (3.20), we also obtain 



eC 



Now, since f" ^ f w-L^(0, 1) we have 



dwo 



dxi 



^0, as e -> 0. (4.6) 



L°°(0,1) 



/"(xi, X2) dx2dxi 



wo(a;i)/'(2;i) dxi 

~^ 9 Wo{xi)fo{xi) dxi, as e — )■ 0. 
Jo 



Therefore, we get from (4.5), (4.6), (4.7) and (2.5) that 



Next we show that 

a^{wo + K*^, Wo + k") — )■ aoiwQ, Wq) as e — )■ 0. 
First we compute the hmit of 

a^iwQ + K^, wq) = / { V(wo + 1^") ■ Vwq + {wq + K'')wo}dxidx2 



(4.7) 



(4.^ 



(4.9) 



-1 



dwp _ Q ^ m X2\ dwp 
dxi ^ V e ' e / dxi 



eX 



xi X2\ d'^wo ^ dwt 



e e / dxf J dxi 



-dxidx2 



+ e / {wq + k'')wq dxidx2 



-1 



dwn"^ 



dxi 



1-dyX 



— e 



eX 



Xi X2 

e ' e . 
xi X2 \ dwo d'^wo 



^ Jr^ 



e e / dxi dx\ 



+ eX 



xi X2 \ dwo 



e e 



Wo 



dxi 



dxidx2 



(4.10) 



as e — )■ 0. Since 



dwo' 



e dxi 



l-dyX 



Xi X2 



e e 



and 



dxidx2 



rg{y) 

/ {l-dyX{y,z)] dz 
Jo 



-^j-^ ~ 9yX ^— , I dzdxi 




JO 
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is a L-periodic function, we obtain that 

e 



^ /o ^ dxi 



Notice that we also have 



WQ\^dxidx2^g I Iwoprfxi as e — 0. 



(4.11) 



(4.12) 



Since wq does not depend on X2, it follows from (3.20) that 
'Xi X2\ dwo d'^WQ 



— e 



+ eX [—, — ) Wo— > 



Hence, we have from ( [4ll| ), ( |4l2| , ( [4l3| and (g that 

«£('?^o + "W^o) 0-0(^0) ""^0)5 as e — 0. 



"^""^ "U^o-T"" dxidx2 — )■ 0, as e — > 0. (4.13) 
e e / ctxi 



(4.14) 



Finally, arguing as in (4.10), we can obtain from (3.20) and (3.21) that 



Oe(tfo + ^^) 0, as e — )• 0, 



getting the statement (4.9). 



Therefore, in accordance with (4.4), we obtain 



completing the proof. 



□ 



Remark 4.2. // the original non homogeneous term does not depend on X2 and e, ie. 



f'^{xi,X2) = f{xi), then it follows from the above definitions (4.1) and (4.3) that 

rtg{xi/e) 



fix, 



f{xi)dx2 = g{xi/e)f{xi) 



w-L\0,l). 



Hence, equation (3.12) is in agreement with the one found via the method of Multiple Scales 
in Section^ 



Remark 4.3. In Remark 3.6 we show the positiveness of the constant r. Hence, the solution 
Wq of the homogenized equation actually exists, is unique and satisfies Wq G iJ^(0, 1) fl 
C\0,1). 



5. Second-order corrector 



Let Wq be the homogenized solution (3.12), X and 6 be the auxihary solutions given by 
(3.7) and (3.16) on the basic cell Y*, which were conveniently defined in the thin domain 
by way off Remark 3.4 
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In this section, we use the second-order corrector (3.19) to present an error estimate when 



we replace the solutions w'^ of (1.2) by the first-order truncation 



with respect to the norm 
second-order truncation 



(5.1) 

\m{R^) introduced in Section [2| Thus, let us consider the 

'Xi X2\ dWo 2 /I /^i d'^WQ 



(5.2) 



Theorem 5.1. Let he the thin domain defined in (1.1) and let he the solution of 
prohlem (fOl) with f'{xi,X2) = f{xi), f G W^'°°{0,1). 



Then, ifWl and are given hy ( |5.1 ) and (5.2) respectively, we ohtain that 

\\\w'-m\\\HHRn<KiV~e. (5.3) 

Consequently, we ohtain the following rate for the first-order approximation 

\\\w'-m\\\HHRn<K2V~e, (5.4) 

where Ki and K2 are positive constants independent of e> 0. 



Proof First, we note that (5.4) is a direct consequence from (5.3) and (3.21), since that 



1 


^2q 


2 




d'^wo 






dx"^ 






dx"^ 


L°°(0,1) 



lfll|2 



+ e 



d'^WQ 



dx'^ 



L°=(0,1) 



for some constant K2 > independent of e. 

Now we estimate the norm 1 1 | = Oe(0^0'') of the function 0^ given by 

In order to do it, we compute 0^(0^, f) for arbitrary test functions (f G H^{R'^) and establish 
an estimate of the form 



\a,{(t)%^)\ < K{e)\\\^\\\m^R^). 
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Since is an elliptic form, we can take = 0^ in the above inequality to obtain the desired 



estimate. Using the notation from Section^ we get by (3.3) that 



i=l 



i=l 



dxx X 



(Pwq 

dx"^ 

dwo 
dx 



+ 2d, 



dwo 
dx 



dx"^ 



Hence, due to (3.7) and (3.16), we have 



r)2 



1=1 



'dx7 



i=l 



+e 



^'"°+2 9„X$| + (l-r-29,X) 



d^ 



Wo 



dx^ 



dx^ 



{X -2dy 



dx^ 



dx"^ 



2 ^ d^wp 
dx^ ' 



Consequently, it follows from (1.2) and ( 3.12[ ) (after some calculations) that 



i=l 



dx-' 



+ 



f + r 



d'^vjQ 




\d^ 


Wo 


dx'^ " 


— e 


dx^ 


« ft'? - 


d'^WQ\ 






dx^ 







{X - 2dy9) - X 



dwQ 
dx 



That is, 



-E 

i=l 



dxi 



+ 



— e 



dx^ 



dx 



— e 



d'^wo d'^wo ' 



dx'^ 



dx^ 



in for all e > 0. 
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On the boundary dR^, we have by the identity (3.4) and boundary conditions from (3.7) 
and (3.16) that 



dm 



dm 



Xl,X2 



wq — eX 



' e 



dx 



dx^ 



- (dr, + -dy] (wo-eX 



dwn n ^ dWn 



dx 



dx"^ 



dJwo-eX^ + eH^\ m 



1 

e 

dwo ( dX 



- 



dx \dN 
dx-^ 



dx 



NA+e 



dwQ 
dx'^ 



dx^ 
XNi 



86 
ON 



m 



dx-^ 



Thus, the function 0*^ satisfies the following boundary value problem 

-A(j)' + = eF^ in R" 

H' on dR' 



dm 



where 

F'{XI,X2) 

and 



d^ 



Wo, 



— e 



dx^ 

Xi X2 





X2^ 


) - 2dye 






/ d^WQ 




d'^wo 


\ dx^ 




dx^ ^ 



Xi X2 



X 



Xi X2\ dwo 
dx 



(xi) 



for a.e. (xi, X2) G R" 



H%x,,x2) = -e(^ 



Xi X2\ dwQ 
dx^ 



Xi) for a.e. {xi,X2) G dR'' 



(5.5) 



(5.6) 



(5.7) 



We consider now the variational formulation of problem (5.5): find 0*^ G H^{R'^) such that 

a,{(j)%ip)= [ F'ipdx + e [ H'N^ipdS. 



(5.8) 



Observe that the function 0*^ must satisfy an uniform a priori estimate on e. Indeed, if we 
take (y9 = 0Mn the expression (5.8), we obtain 



wmwmiRn 



(5.9) 



- |ae(0^0^)| 

< ||01U2(k<!)||-^^||l2(/?o + e \\4>^\\L^{dRA\H^Nl\\L2(dR^). 

We need to get sharp inequalities on F*^ and H"^ to estimate 0^(0*^,0^). It is clear from 
their definitions that these estimates will be consequence of those ones for wq, X and 9. 
Since / is a smooth function, we have by classical regularity results given in [1] that the 
solution Wo of the homogenized problem is smooth enough to guarantee that its derivatives 
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up to the fourth order are in L°°(0, 1). Note that similar statements are also true for X and 
9 e H\Y*). 

Due to the periodicity of X, we have by (3.20) that 
IIX 



.^ < X 



< 



)Jj^\\0\\l^y*) and | |9y6'| |l2(r.) < ^J^\\6\\l^y*)- 



Consequently, it is clear from (5.6) that there exists Kq independent of e such that 



(5.10) 



Let us observe that Kq depends on the period L of the norms of X, 9 and dyO in L^(y*), as 
well of the norms of ^t^, and ^ in L°°(0, 1). 

Now, let us denote the oscillatory part of OR" by doR^ = {{xi,eg{xi/e)) : < xi < 1}, the 
fixed part by dfR'' = {{xi,0)) : < a;i < 1} and the lateral part of dW as diW = {(0,a;2)) : 
< X2 < ^g{0)} U {(1,X2)) : < a;2 < eg{l/e)}. From definition (5.7) we have 



9 



< 



< 



< 



< 



dx^ 



Xl X2\ dwn 



e ' e / dx"^ 



-{xi) Nl{xi,X2) 



dS 



L°°(0,1) JdR' 



dwQ 



dx^ 



L°°{0,1) \JdoR' 



Xl X2 

e e 

Xl X2 

e e 



dS 



dS + 



dfR'^ 



Xl X2 



e e 



dS 



\fc=i k=i 



~j^\r\\L^{aY*) 

,3 /^^3| 



where Ki = ||(1 + g')^\\L'^(o,L) \\dw'o/dx \\j^^^q -^-^ is independent of e. Note that we have 

used the periodicity of 6 to get J^^^, [^("t", ^)| dS = 0. Consequently there exists K2 > 
independent of e such that 

Now we have all the ingredients to estimate a 



(5.11) 



Due to (5.10) and (5.11) we get 



from (5.9) that 



e| I |2 



||0l|L2(iJ^) +(^Ko ||0'||L2(a^.). (5.12) 

Hence, the desired result follows from the following fact: If (/? G H^lR^), then there exists a 



constant C independent of e such that 



(5.13) 



Indeed, if we combine (5.12) and (5.13), we obtain Ki > independent of e such that 

\\m{R') < Ki e^/^. 
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The proof of (5.13) can be found in [T2[[Td]. We recall it here for the reader's convenience. 
From smoothness of BY* we can define a smooth extension, M = (Mi,M2) G C^{Y*), of 
the unitary normal vector field on Y*, such that M{y,z) = N{y,z) a.e. dY*, and with 
support of M in a some neighborhood of dY*. Hence, for all (f G H^{R'^) it follows that 

= y V (<^2) . M ^) dx,dx2 + j divM ^) dx,dx2 



2 



j ^V^.M(^^,^yxidx2 + e-' j y,^I^J2dyM{^,^)'^dx,dx2 



< Ci\\v\\l2{r-)\\'^v\\l2(r^) + C2e ^W^Wh^R.) 

< max{Ci,C2}e"^||(y9||^i(^,). 



□ 



6. Final conclusion 



In this work, we give a precise rate of the convergence for solutions of an elliptic problem 
posed in a family of rough domains with a singular collapsing structure. In our analysis we 
use the corrector approach grounded in a formal asymptotic expansion of solutions, widely 
used in homogenization theory, to obtain a rigorous strong convergence result in Theorem 



4.1 The second-order corrector is used to obtain an error estimate for the convergence result 



in Theorem 5.1 We observe that the rate of our convergence result is sharp in the sense 



that we do not have introduced any boundary layer terms in the correctors, see e.g. [21 



Acknowledgments. Part of this work was done while the authors were visiting the Depar- 
tamento de Matematica of the Universidade de Sao Paulo, SP - Brazil. We kindly express 
our gratitude to the Department. The authors thanks to Jose M. Arrieta and Alexandre N. 
Carvalho for their remarks and suggestions. We also would like to thank the anonymous 
referee whose comments have considerably improved the writing of the paper. 

References 

[1] S. Agmon, A. Dougiis, L. Nircnberg; Estimates near the boundary for solutions of elliptic partial diffe- 
rential equations satisfying general boundary value conditions I; Comm. Pm'e Appl. Math. 12 (1959) 
623-727. 

[2] Y. Amirat, O. Bodart, U. de Maio, A. Gaudiello; Asymptotic Approximation of the solution of the 
Laplace equation in a domain with highly oscillating boundary; SIAM J. Math. Anal. 35, (2004) 1598- 
1616. 

[3] J. M. Arrieta; Spectral properties of Schrodinger operators under perturbations of the domain; Ph.D. 

Thesis, Georgia Inst, of Tech. (1991). 
[4] J. M. Arrieta, A. N. Carvalho, M. C. Pereira and R. P. Silva; Nonlinear parabolic problems in thin 

domains with a highly oscillatory boundary; Nonhnear Analysis: Theory Methods and Appl., (74) 15 

(2011) 5111-5132. 



ERROR ESTIMATE FOR A NEUMANN PROBLEM 



17 



[5] J. M. Arrieta and M. C. Pereira; Elliptic problems in thin domains with highly oscillating boundaries; 
Bol. Soc. Esp. Mat. Apl. 51 (2010) 17-25. 

[6] J. M. Arrieta and M. C. Pereira; Homogenization in a thin domain with an oscillatory boundary; J. 
Math. Pures et Appl. 96 (2011) 29-57. 

[7] A. Bensoussan, J. L. Lions and G. Papanicolaou; Asymptotic Analysis for Periodic Structures; North- 
Holland (1978). 

[8] R. Brizzi, J. P. Chalot; Boundary homogenization and Neumann boundary problem; Ricerce di Matem- 

atica XLVI, 2 (1997) 341-387. 
[9] D. Cioranescu, A. Damlamian and G. Griso; The periodic unfolding method in homogenization; SIAM 
J. Math. Anal. 40 no. 4 (2008) 1585-1620. 
[10] D. Cioranescu and P. Donato; An Introduction to Homogenization; Oxford lecture series in mathematics 

and its applications (1999). 
[11] D. Cioranescu and J. Saint J. Paulin; Homogenization in open sets with holes; J. Math Anal. Appl. 71 
(1979), 590-607. 

[12] D. Cioranescu and J. Saint J. Paulin; Homogenization of Reticulated Structures; Springer Verlag (1980). 

[13] A. Damlamian, K. Pettersson; Homogenization of oscillating boundaries; Discrete and Continuous Dy- 
namical Systems 23, (2009), 197-219. 

[14] T. Elsken; Continuity of attractors for net-shaped thin domain; Topol. Meth. Nonlinear Analysis, 26 
(2005) 315-354. 

[15] J. K. Hale and G. Raugel; Reaction- diffusion equation on thin domains; J. Math. Pures et Appl. (9) 71 
no. 1 (1992) 33-95. 

[16] J. -L. Lions; Asymptotic expansions in perforated media with a periodic structure; Rocky Mountain J. 

Math. 10 (1) (1998) 125-140. 
[17] D. N. Arnold and A. L. Madureira; Asymptotic Estimates of Hierarchical; Modeling Mathematical 

Modeling and Methods in Applied Sciences, vol. 13, No. 9, (2003) 1325-1350. 
[18] A. L. Madureira and F. Valentin; Asymptotics of the Poisson Problem in domains with curved rough 

boundaries; SIAM Journal on Mathematical Analysis 38, No. 5, (2007) 1450-1473. 
[19] T. A. Mel'nyk; Homogenization of the Poisson equation in a thick periodic junction; Z. Anal. Anwen- 

dungen 18, (4), (1999) 953-975. 
[20] J. Nevard and J. B. Keller; Homogenization of rough boundaries and interfaces; SIAM J. Appl. Math., 

57 no. 6 (1997) 1660-1686 
[21] G. Panasenko; Multi-scale modelling for structures and composites; Springer, Dordrecht, 2005. 
[22] M. Prizzi and K. P. Rybakowski; The effect of dom.ain squeezing upon the dynamics of reaction- diffusion 

equations; Journal of Diff. Equations, 173 (2) (2001) 271-320. 
[23] M. Prizzi and M. Rinaldi and K. P. Rybakowski; Curved thin domains and parabolic equations; Studia 

mathematica, 151 (2) (2002) 109-140. 
[24] G. Raugel; Dynamics of partial differential equations on thin domains; Lecture Notes in Math., vol 

1609, Springer Verlag (1995). 
[25] E. Sanchez-Palencia: Non-Homogeneous Media and Vibration Theory; Lecture Notes in Phys. 127, 

Springer Verlag (1980). 

[26] R. P. Silva; Semicontinuidade inferior de atratores para problemas parabolicos em dommios finos; Phd 

Thesis, ICMC - USP, (2007). 
[27] L. Tartar; The General Theory of Homogenization. A personalized introduction; Lecture Notes of the 

Un. Mat. Ital, 7, Springer- Verlag, Berlin (2009). 



18 



M. C. PEREIRA AND R. P. SILVA 



(Marcone C. Pereira) Escola de Artes, Ciencias e Humanidades, Universidade de Sao Paulo, 
Sao Paulo SP, Brazil 

E-mail address: marcone@usp.br 

(Ricardo P. Silva) Instituto DE Geociencias e Ciencias Exatas, Universidade Estadual Paulista, 
Rio Claro SP, Brazil 

E-mail address: rpsilva@rc.unesp.br 



